Prospects for direct detection of primordial gravitational waves. 
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We study the primordial gravitational wave background produced in models of single field in- 
flation. Using the inflationary flow approach, we investigate the amplitude of gravitational wave 
spectrum, Wgw, in the frequency range 1 mHz - 1 Hz pertinent to future space-based laser interfer- 
ometers. For models that satisfy the current observational constraint on the tensor-to-scalar ratio, 
r < 0.36, we derive a strict upper bound of Ug-^j < 1.6 x 10~^^ independent of the form of the 
inflationary potential. Applying, in addition, the observational constraints on the spectral index ria 
and its running, Wgw is expected to be considerably lower than this bound unless the shape of the 
potential is flnely tuned. We contrast our numerical results with those based on simple power-law 
extrapolation of the tensor power spectrum from CMB scales. In addition to single fleld inflation, 
we summarise a number of other possible cosmological sources of primordial gravitational waves and 
assess what might be learnt from direct detection experiments such as LISA, Big Bang Observer 
and beyond. 



I. INTRODUCTION 

The existence of a stochastic background of primordial 
gravitational wave from inflation has yet to be verified by 
observation. A significant detection would not only con- 
firm the success of inflation, but would also serve as a 
unique observational window to physics during the very 
early universe. Since the first resonant bar of Joseph We- 
ber in the 1960s ll|, direct-detection experiments such as 
LIGO have reached the stage where detection of astro- 
physical sources is a realistic prospect. Discussion of am- 
bitious space-based interferometers beyond LISA is well 
underway (see Table I for summary and references). One 
of the main goals of post-LISA missions is to detect the 
stochastic gravitational wave background predicted by 
inflation. The most ambitious of these proposed exper- 
iments looks forward to a precision limited only by the 
Heisenberg uncertainty. 

In the context of inflationary models, the amplitude of 
the stochastic gravitational wave background remains ex- 
tremely uncertain because neither the energy scale of in- 
flation, nor the shape of the inflaton potential, is known. 
Previous studies ^, ^ have often relied on some form of 
potential to calculate the gravitational wave spectrum. 
While a fuller understanding of the inflationary mech- 
anism (if indeed inflation occurred) awaits further de- 
velopment in fundamental physics, we ask what generic 
predictions, relevant to direct gravitational wave experi- 
ments, can be made in simple models of inflation without 
recourse to speciflc potentials. In this paper, we address 
this problem and assess the future prospects for direct de- 
tection experiments as they confront inflation and other 
theoretical ideas. 

After a brief overview of inflation, we calculate the 
amplitude of primordial gravitational wave spectrum pre- 



dicted by inflation and comment on the main uncertain- 
ties involved in this calculation. We then generate models 
of inflation stochastically using the inflationary flow ap- 
proach and study the gravitational wave amplitudes in 
these models. Direct detection experiments probe physi- 
cal scales that are at least 15 orders of magnitude smaller 
than the scales probed by CMB experiments. The infla- 
tionary flow approach allows us to investigate the limita- 
tions of simple extrapolation between these scales using 
a 'slow-roll' approximation. Next, we briefly discuss a 
range of other mechanisms, in addition to single field in- 
flation, for generating primordial gravitational waves at 
direct detection scales. Finally, we assess the prospects 
that future gravitational wave experiments might shed 
light on inflation and the early universe. 



II. INFLATIONARY PERTURBATIONS 

We shall work in the so-called "Hamilton-Jacobi" for- 
mulation, in which the Hubble parameter H describes 
the inflationary dynamics. The 'slow-roll' parameters e 
and T] are defined in terms of the inflaton-valued Hubble 
parameter i?(0) as follows: 
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where primes denote derivatives with respect to the in- 
flaton value (f) and mpi is the Planck mass. Following 
the normalizations of the amplitudes A{k) of pri- 

mordial power spectra V{k) are given, to lowest order, 
by 
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where S and T denote scalar and tensor components 
respectively. The amplitudes are evaluated when each 
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Experiment 


Time-scale 


Sensitivity to ojgw 


Optimum Frequency (Hz) 


Reference 


Advanced LIGO 


2009 


10-8 


100 


\£ 


LISA 


2014 


10-" 


0.005 




BBO/DECIGO 


2025? 


10"" - 10-^^ 


0.1 




Ultimate DECIGO 


2035? 


10-2" 


0.1-1 





TABLE I: Summary of some relevant parameters of future experiments for direct detection of gravitational waves (here 
iOgw = flgvihg). These experiments include the ground-based Advanced Laser Interferometer Gravitational- Wave Observatory 
(LIGO), as well as space missions such as the Laser Interferometer Space Antenna (LISA), NASA's Big Bang Observer (BBO) 
and Japan's Deci-Hertz Interferometer Gravitational Wave Observatory (DECIGO). The ultimate DECIGO is envisaged to be 
a quantum limited interferometer in space with 100-kg test masses. The quoted time-scales and sensitivities are indicative only. 



mode, fc, is equal in scale to the Hubble radius, i.e. when 
k = aH. As the inflaton evolves, the rate at which differ- 
ent scales leave the Hubble radius is given by 
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Small departures of the primordial spectra from scale 
invariance are measured by the spectral indices defined 
as 
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In practice, however, it is common to let the spectral 
indices quantify variations around a pivot scale kg. In 
this approximation, the power spectra are parametrized 
by: 



Vsik) 



Vrik) = Vriko) 



ko 



(7) 



(8) 



Using Equation one finds that the spectral indices 
can be approximated to 0{e,r]) by 

(9) 
(10) 



ris - 1 ~ 2ry - 4e 
riT — — 2e . 



Often, it is convenient to describe a power spectrum as 
blue when its index exceeds unity, or red otherwise. In 
this terminology, the tensor power spectrum is said to 
always be tilted red. However, Pre-Big Bang and cyclic 
scenarios provide exceptions, where the tensor spectrum 
is strongly blue. We return to this point in Section Ivl 

The ratio between the tensor and scalar amplitudes is 
clearly 



In concordance with Ref. 
tensor-to-scalar ratio r as: 
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Equations H1(J|) and (|12|l combine to give the lowest order 
consistency relation: 
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Note that the definition of r varies widely in the lit- 
erature. For instance, it is often defined as the ratio 
of tensor to scalar quadrupole CMB anisotropy r2 = 
(|a2mP)/(|a2mP) QEI - Such a definition is cosmology- 
dependent, especially on the dark energy density f^A- 
The conversion is lla: 



r2 



0.84- 0.025f7A 



1.04-0.82fiA 



2f7i 



(14) 



III. GRAVITATIONAL WAVE SPECTRUM 

We now briefiy derive an expression for the primor- 
dial gravitational wave spectrum in terms of infiationary 
observables r, Us and nx- This Section establishes the 
definitions and normalizations of various quantities used 
in the rest of the paper. This is important because there 
are a number of derivations of the gravitational wave en- 
ergy spectrum expected from inflation in the literature, 
of varying accuracy. The discussion here is based on Refs. 

iniiiiiiiia 

We begin by considering the primordial gravitational 
waves produced via tensor perturbation hij of the 
flat Friedmann-Robertson- Walker metric. In the syn- 
chronous gauge (ft,po = 0) and natural units {c — h—1), 
the perturbed metric is 



ds^ = ~dt^ -I- a^(t){Sij + hij)dx'dx^ 



(15) 



where a{t) is the scale factor in coordinate time. By 
further imposing the transverse traceless conditions, the 
tensor perturbations can be described by two polariza- 
tion states h\{x,t) with A = +, x. In Fourier space, the 
tensor power spectrum Pxik) observed today (t — to) is 
given by the variance 
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Relative to the background FRW cosmology, an effec- 
tive stress-energy tensor of g ravitational waves can be 
defined unambiguously as [2l| 
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The component — Tq = gives the energy density 
of gravitational wave background. 
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The strength of the primordial gravitational waves is 
characterized by the gravitational wave energy spectrum: 
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where pc — 3Hq/8ttG is the critical density and Hq = 
lOO/io kms"^Mpc"^ Substituting into (HHJ gives an im- 
portant result: 
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which is consistent with Ref. [23|. The physical density 
in gravitational waves is defined as 



(21) 



and is independent of the value of Hq. Following previous 
work we shall calculate constraints on the quantity 

Next, ignoring anisotropic stresses, the Einstein equa- 
tions require that each state hx{k) evolves via the mass- 
less Klein- Gordon equation 



d^hx 2 da dhx 
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where r is the conformal time. Anisotropic stresses from 
free streaming particles can create a non-zero source term 
on the right-hand side of Equation H22I) . We return to 
this point shortly. 

The tensor power spectrum at the end of inflation, 
Prik), can be related to the tensor power spectrum at 
the present day by a transfer function T{k) , 



PT{k)^T\k)PT{k) 
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By numerically integrating Equation l|22(l . the transfer 
function is found to be well approximated by the form 
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where kcq = O-Ol^flmh^ Mpc~^ is the wavenumber cor- 
responding to the Hubble radius at the time that mat- 
ter and radiation have equal energy densities. Using the 
cosmological parameters determined by combining data 



from several surveys [2j|, one finds k^q = 0.0104 Mpc"-*^ 
and To = 1.41 X 10'* Mpc. Combining Equations 
and H24|l gives the gravitational wave energy spectrum 
for k ^ kc 



r2gw(fc) 
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At present, the best constraints on the normalization 
of the tensor spectrum come from CMB anisotropy ex- 
periments. It is tempting therefore to evaluate Equation 
(|25|l by normalizing at CMB scales. However, the physi- 
cal scales probed by CMB experiments are about 15 or- 
ders of magnitude larger than the scales probed by direct 
gravitational wave detection experiments. In the context 
of this paper, there are both positive and negative aspects 
associated with this large difference in scales. On the one 
hand, it is not straightforward to extrapolate from CMB 
scales and infer what might be observed by direct detec- 
tion experiments, even under the restrictive assumption 
of single field inflation. On the other hand, this large dif- 
ference in scales means that direct detection experiments 
offer the prospect of learning something fundamentally 
new that cannot be probed by CMB experiments. The 
main aim of this paper is to investigate how reliably one 
can extrapolate Equation (|25|l from CMB scales, with as 
few constraints on the form of the inflationary potential 
as possible. 

Although a tensor component has not yet been ob- 
served in the CMB anisotropics, the amplitude of the 
scalar component has been determined quite accurately. 
At a fiducial 'pivot scale', ko = 0.002 Mpc~^, the com- 
bined results from WMAP, 2dFGRS and Lyman a sur- 
veys give |22 



Ps{ko = 0.002 Mpc"*) ~ 2.21 x 10" 
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Using the above result and expressing k in terms of phys- 
ical frequency / = fc/27r, we finally obtain an expression 
for primordial gravitational wave spectral energy in terms 
of / and inflationary observables r and nx only: 
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where fo = 3.10 x 10~^^ Hz. This relation is valid as long 
as / ^ /cq ^ 10^*'' Hz and ut is independent of scale. 

Further, if r and nr are accurately approximated by 
first order expressions in e. Equation H27() becomes 



Wgw(/) 6.98 X IQ-^^el^ 
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This expression is maximized at e = [21n(///o)] *, with 

6.98 X lO-i-* 
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According to this approximation, the strength of primor- 
dial gravitational waves at direct detection scales does 
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not increase proportionally with r because models with 
large r have a large red tensor tilt. The crucial assump- 
tion is, of course, that the power-law parametrization 
Vxik) oc fc"^, with constant index ht, remains accurate 
over the many orders of magnitude from CMB scales to 
those probed by direct detection experiments. In the 
next Section, we use numerical calculations of inflation- 
ary evolution to go beyond this approximation, finding 
many examples of inflationary potentials for which Equa- 
tion 1)27(1 is violated badly. 

Finally, we comment on suggestions that tensor power 
may be significantly reduced due to anisotropic stresses 
from free-streaming neutrinos [23L . For three stan- 
dard species of neutrinos, Wgw is damped by a factor 
of < (0.80)^ on scales which re-entered the Hubble ra- 
dius during radiation era after neutrino decoupling at a 
temperature of a few MeV. These scales correspond to 
frequencies of about 10~^^ Hz, well below the frequen- 
cies relevant to direct detection of gravitational waves. 
Damping at direct detection frequencies is still possi- 
ble via more complicated mechanisms, for instance, free- 
streaming of exotic massive particles which decouple 
above the electroweak scale, or perhaps via extra dimen- 
sional physics manifesting above the TeV scale (see Sec- 
tion 01. But because these phenomena are still specu- 
lative and poorly understood, we have chosen to ignore 
them at present. For a review of these and other damping 
mechanisms, see Hfl- 



IV. NUMERICAL METHOD 

As we have discussed above, it is interesting to anal- 
yse the stochastic gravitational wave background with- 
out relying on specific forms for the inflaton potential. 
Given our lack of knowledge of the fundamental physics 
underlying inflation, we have tackled this problem by in- 
vestigating a large number of viable inflationary models 
numerically. 

Our approach is based on the inflationary flow equa- 
tions, first introduced in Ref [2^ and further developed 
in [H m m mm in the notation of Ref. 



the flow equations are: 
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Here the derivative with respect to the number of e-folds, 
N, runs in the opposite direction to time. The flow equa- 
tions represent an infinite dimensional dynamical system 
whose dynamics is well understood jlS]. The parame- 
ters of the system are given in terms of inflaton-valued 



Hubble parameter H{(j)) by: 

An J 
2?7 - 4e . 



H 



_ mpi 
47r 

-1 d'+^H 



El 
IT 



(31) 



The hierarchy completely defines the function H{(j)), 
which in turn determines the inflaton potential V{cj)) via 
the Hamilton-Jacobi Equation, 
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In terms of the flow parameters, the inflationary observ- 
ables are given to next to leading order by j^S] 
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where C = 4(ln 2 + 7) - 5 ~ 0.0814514 (with 7 the Euler- 
Mascheroni constant). Variations of the spectral indices 
with scales are approximated to first order by the "run- 
nings" dns/dlnk and dnx / dink. While ut may be mea- 
sured directly by BBO/DECIGO (via the slope of Wgw 
around 1 Hz |3l|) or indirectly (via the consistency re- 



lation j32l l33|). its running, however, is likely to remain 
poorly constrained in the foreseeable future. Thus, we 
have not explicitly analysed the gravitational wave spec- 
trum with respect to dnx/dhik. 

We ran a program (previously used in 0, 0]) that 
generates models of inflation stochastically. The program 
flrst selects the initial conflguration of a model from uni- 
form distributions within the following ranges: 

eo e [0, 0.8] , 
<Jo e [-0.5,0.5] , 

Co e [-0.05,0.05] , (36) 
^Xh\o e [-0.025 X 5"^+^ 0.025 X 5-^+^] , (3 < ^ < 10) 
"AhIo = , 

where the hierarchy is truncated at £ = 10. Each model 
is evolved forward in time (backward in e-fold) until in- 
flation ends in one of the following ways: 

1. By achieving e — 1. When this happens, we say for 
convenience that the 'slow-roll' condition has been 
violated. Observables on CMB scale are then cal- 
culated 60 e-folds before the end of inflation. This 
number of e-fold at which observables are generated 
is in accordance with the analyses of Refs. |3i US ■ 
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FIG. 1: (Colour online) Plots of gravitational wave spectrum ujgw against tensor-to-scalar ratio r for a large number of models 
evolved with the inflationary flow equations. Square (red) points indicate models satisfying the observational constraints on 
ris and dus/dlnk given by 1431 . In panel (a), ojgw is calculated using the extrapolation formula 12711 . The solid line is the first 
order approximation given by Equation 1281 . In panel (b), Wgw is calculated using the formula 14111 . The solid (green) curve in 
panel (b) shows the bound given by Equation 14211 . with the parameter A = 7 . 



2. By an abrupt termination, perhaps from interven- 
tion of an auxiliary field as in hybrid inflation [36| | , 
or, when open strings become tachyonic in brane in- 
flation js^, m, 113 • Because these scenarios accom- 
modate a large number of e-folds during inflation, 
one identifies them with an asymptotic behaviour of 
a trajectory. In practice, those models inflating for 
more than 200 e-folds are grouped under this cat- 
egory. The observables are then calculated along 
the asymptote. 



We 
model 



produced 
calculated 



10^ realizations and for each 
five key observables, namely 



{r, Us, riT, dus/dlnk, cugw}- Working with next 
to leading order expressions in e,ry, we use the following 
expression for the primordial tensor power spectra with 
the assumption that e and rj are approximately constant 
as each mode crosses the Hubble radius jlQj 
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where C is defined as before. The gravitational wave 
spectrum depends on Equation evaluated when the 
direct detection scales cross the Hubble radius. Since 
modes with frequencies in the direct detection range of 
around 0.1-1 Hz exit the Hubble radius when TV ~ 20, 
the relation between the Hubble parameters at direct de- 



tection and CMB scales is given by 

HcMBexp(^-J^ e{N)dN^ . (38) 

The gravitational wave spectrum is now given in terms 
of the flow-parameters at scale ko by: 

15 



Vs{ko)rIik) 



(39) 



with r given by Equation H33|) and 
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Inserting numerical factors gives: 

Wgw(fc) 



4.36 X IQ-^^r. 
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For comparison between Equation H41|) and the extrap- 
olation formula H27|l . we evaluate the gravitational wave 
spectrum in our models using both expressions. We 
adopted a nominal BBO/DECIGO frequency of 0.1 Hz, 
consistent with Ref. In any case, the results are in- 
sensitive to the choice of frequency as long as the latter 
exceeds the neutrino damping scale {^-^ 10~^^ Hz). 



A. Dependence of cjgw on r 

Figure ^ summarizes our main results. Most of the 
models are of the 'hybrid' type for which the tensor mode 
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FIG. 2: (Colour online) Some trajectories H{N), from CMB scales (A'^ ~ 60) to the end of inflation (A'^ — 0), for models 
evolved using the inflationary flow equations. The models plotted in panel (a) have high gravitational wave amplitudes at direct 
detection scales (ugw > 2.5 x 10~^®), whilst those shown in panel (b) have low amplitudes (ugw < 5 x 10~^^). All of these 
models satisfy the observational constraints on Us and dris/dlnfc given by Equation 1431 . and have high tensor amplitudes in 
the range 0.15 < r < 0.25. 



is negligible (r w 0, Wgw ~ 0) and in which the stochastic 
gravitational wave background is well below the detec- 
tion threshold of any conceivable experiment. Figure^ 



shows the results in the 



plane when the extrap- 



olation formula (|27|1 is used to compute c^gw Most of 
the 'non-trivial' models {i.e. models with high Wgw) lie a 
few percent below the first order prediction (|28|l shown 
by the solid (green) line. All of these non-trivial models 
achieve e = 1 at the end of inflation. Fig. QJ) shows the 
results of using the formula (|41|l to compute Wgw. The 
distribution of points now spans a large fraction of the 
r — Wgw plane. The inflationary flow formulation shows 
that the first order extrapolation formula H27|) is too re- 
strictive. Since the shape of the inflationary potential is 
unknown, it is not possible to extrapolate reliably from 
CMB scales to the much smaller scales probed by direct 
detection experiments. Figure ^ shows that it is possi- 
ble to find inflationary models in which, for instance, the 
flow variables change rapidly within the last e- folds, thus 
enhancing ujg^ at direct detection scales. 

The solid (green) line in Figure^ shows the expres- 
sion, 



4.36 X 10"^^r 



1 - 
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(42) 



where the constant A ~ min(e(A:o)/e(A:)) depends on the 
distribution In our runs, we find A ^ 7. This ex- 

pression provides an accurate upper bound to Wgw. Equa- 
tion H42|) simply expresses the constraint that the Hubble 



parameter is constant between CMB and direct detection 
scales, modulated by the term in square brackets which 
expresses the details of how inflation ends. However, for 
any value r < 1, the term in square brackets is close to 
unity and so is insensitive to the parameter A and hence 
to the distribution Ij^til) . 

The (red) square points in Figure ^ show the subset 
of models that satisfy the 2a observational constraints 
[23, oil and dris/dlnk, 



0.92 < Us < 1.06, 



-1.04 < 



dink 



< 0.03. 



(43) 



These models roughly follow the locus of the first order 
extrapolation shown in Figure^, but with a large scat- 
ter. As a conservative bound we apply Equation H42() 
with the observational constraint r < 0.36 [23, to give 



We 



< 1.6 X 10 
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(44) 



As this paper was nearing completion, a paper by p4| 
appeared describing a similar analysis. Our results are 
broadly compatible, but there appear to be some discrep- 
ancies. Comparing our Figure ^ with their Figure 2, 
we see that the swathe of points satisfying H43|) matches 
roughly the shape of the contoured region in their Figure. 
However, we find models with low values of Wgw < 10"^^ 
at all values of r whereas they do not. Furthermore, at 
high values of r > 0.1, they appear to find models that 
lie above the bound given by H42|) . Their results do not 
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FIG. 3: (Colour online) The gravitational wave spectrum cugw plotted against scalar spectral index for a large number of 
models evolved using the inflationary flow equations. Square (red) points indicate models satisfying the observational constraints 
on Us and dna/dlnk (Equation 14311 1 and satisfying r < 0.36. In panel (a), Ljgw is calculated using the extrapolation formula 
1271 1. The (green) solid curve shows the bound given by Eq. 1451 1. In panel (b), tjgw is calculated using formula 14H and the 
flow equation integration. 



seem physically plausible to us ^ . 

Examples of some trajectories H{N), from CMB scales 
to the end of inflation, are shown in Figure |21 AH of 
these models satisfy the observational constraints on Ug 
and dus/dhik of Equation H43|l and, in addition, we have 
imposed the constraint 0.15 < r < 0.25, i.e. the mod- 
els have high tensor amplitudes. The models plotted in 
Figure have high gravitational wave amplitudes at 
direct detection scales {ujg^ > 2.5 x 10"^^). In these 
cases, the Hubble parameter stays almost constant from 
= 60 to iV = 20 but declines rapidly thereafter. In 
contrast, the models shown in Figure l^b have low am- 
plitudes ujgvj < 5 X 10"^''. In these cases, H{N) de- 
clines more rapidly between N — 60 and N = 20. These 
sample trajectories show that models with sharp features 
in H{N) (and hence also in V{(j))) within the last 20 
e-folds of inflation will be the first to be ruled out by 
BBO/DECIGO-type detectors. 



B. Dependence of ci;gw on scalar tilt Us 



Figure O shows the models plotted in the Ug 
plane. The extrapolation method (Fig. places most 



of the 'non-trivial' models within a vertical band centered 
around ~ 0.8. The band is sharply capped by the solid 
(green) curve given by differentiating Eq. : 



= 4.24 X 10" 



17.235- 1.303n, 



2.565- 0.541ns 



(45) 



FigEt* shows the distribution when the flow formula- 
tion (|4I|I is used to calculate cjgw. The region beyond the 
envelope (|45l) is now populated by many models, some of 
which produce Wgw in excess of 10"^''. However, all of 
the models with such high values of are inconsistent 
with the observational constraints on and dris/dhik. 
The (red) squares in FigureOlindicate models that satisfy 
the 2ct observational constraints of Equation (|43|l . and, 
in addition, have r < 0.36. The vast majority of these 
models lie below the line defined by Equation 145() . How- 
ever, it is possible, though rare, for models satisfying the 
observational constraints to exceed Wgw > 3 x 10^^^, 
as given by Equation H45|l . Evidently, one can see from 
Figure^^ that no model satisfying the observational con- 
straints can exceed our conservative bound H44|) . 



V. PROSPECTS FOR DIRECT DETECTION 



^ This is because Ref. |41| extrapolates from CMB scales to direct 
detection scales (using nx and its running) in order to test the 
consistency relation. We thank Hiranya Peiris for clarification. 



The results of the preceding Section show that simple 
single-field inflation models must satisfy the conserva- 
tive constraint of Equation (|44|l at direct detection scales. 
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Furthermore, unless the inflationary parameters are spe- 
ciahy tuned, most single- field inflation models will pro- 
duce Wgw < 3x 10-16. Thus, at the BBO/DECIGO sensi- 
tivities of ^ 10^1^ — 10^^ (see Table 1), a direct detection 
of a stochastic background of gravitational waves would 
be expected only if the inflationary potential contains a 
feature at '--^ 20, as shown in the trajectories plotted in 
Figure 13 This is true even if the tensor-to-scalar ratio is 
high at CMB scales. This is the main conclusion of this 
paper. 

Although this may seem a somewhat pessimistic con- 
clusion for direct detection experiments, it is worth men- 
tioning a range of other cosmological sources (summa- 
rized in Table II) that could produce a stochastic back- 
ground of gravitational waves at direct detection scales. 

(i) Pathological potential: A sudden decrease in energy 
scale of the inflationary universe could be attributed to 
a first order phase transition brought about by the spon- 
taneous symmetry breaking of a field coupled to the in- 
fiaton. As a result, the potential also acquires a 

sharp feature in the form of steps 0jE3j0.EMj kinks 
[4^ or combination of these at various scales . 

In particular, the primordial gravitational wave ampli- 
tude in the so-called 'broken-scale invariance' models has 
been considered in Refs. [^l50l|. which found roughly an 
order of magnitude increase above that given by Equation 
(|28|l . Clearly, a first order phase transition has a negli- 
gible enhancement effect on modes at direct detection 
scales unless the transition occurs at late stages (within 
the last ~20 e-folds) of inflation. On the other hand, if 
scale invariance is broken at around the CMB/LSS scales, 
as suggested by |^ Is^, , then the gravitational wave 
amplitude may be enhanced at scales probed by the fu- 
ture CMB polarization experiments. 

(a) Bubble nucleation: A phase transition may also be 
accompanied by a rapid nucleation of vacuum bubbles 
|54l I55I Is^ . which upon collision during inflation pro- 
duce a large gravitational wave background with of 
order ~ 10"'^ around the direct detection frequencies. 
However, bubble collision at a much lower energy, e.g. 
the electroweak scale, produces virtually negligible grav- 
itational waves with ojgw of order lO"^"^ 57]. In super- 
symmetric extentions of the standard model, this value 
may be larger by several orders of magnitude [s^ and 
perhaps as large as ^ lO^^^ for some parameter choices 
in next-to-minimal models. 

(Hi) Turbulence: A large injection of energy into the cos- 
mological plasma following bubble collision could also set 
up a Kolmogorov spectrum of turbulence. Calculations 
in Refs. |57l Is^ estimate the gravitational wave back- 
ground from turbulence to be comparable to that from 
bubble nucleation. If the turbulence is sourced also by a 
helical field {e.g. primordial magnetic fields), a secondary 
contribution of Wgw ^ 10"^^ is predicted at direct detec- 
tion scales . Relation between Wgw and the strength of 
primordial magnetic fields is further discussed in (61i I63 | 



(iv ) Cosmic strings: A stochastic network of strings 
|63l Ig^ produces a gravitational wave spectrum with a 
long plateau extending from / ^ 10^^° Hz across di- 
rect detection scales [63, Although CMB observa- 
tions show that strings cannot be solely responsible for 
structure formation 67], they can arise in certain models 
of hybrid and brane infiation as a sub-dominant contri- 
bution to the fiuctuations 38, 68, 6_3j]. Recently Refs. 

have calculated the gravitational wave spectrum 
from bursts associated with cusps and kinks in loops of 
cosmic (super) strings as a function of the theoretically 
uncertain intercommutation probability. They conclude 
that the gravitational wave bursts from strings with ten- 
sions as low as G^i ^ IQ^^^ could result in Wgw as large 
as ^ 10~ii. This is potentially detectable by LISA and 
may even be observable by LIGO if Gji > 10"^° and the 
intercommutation probablity small. 

While inflation may be accompanied by all of the phe- 
nomena mentioned above, some alternatives to slow-roll 
inflation have altogether different predictions regarding 
the production of primordial gravitational waves at di- 
rect detection scales. 

(v) Pre- Big Bang and cyclic models: In Pre-Big Bang 
scenarios [t^, E3i a dilaton-driven phase with H > 
gives rise to a gravitational wave amplitude which in- 
creases with frequency f^) for all modes exiting the 
Hubble radius during the Pre-Big Bang era. The pri- 
mordial tensor spectrum in this case is strongly blue 
with ut = 3. The gravitational wave spectrum could 
peak at direct detection scales with amplitude Wgw as 
high as 10~6. within reach of advanced terrestrial detec- 
tors mizElzil- When combined with CMB polariza- 
tion experiments, a strongly blue tensor spectrum can 
be easily ruled out. Nevertheless, the prediction of such 
a large gravitational wave amplitude at direct detection 
scales is sensitive to physics during the 'bounce' around 
t = 0, which remains poorly understood Wt\. In contrast, 
the cyclic model 78, 79] predicts a blue tensor spectrum 
(riT = 2) but with negligible grav itational wave ampli- 
tude at direct detection scales [S^l- 

(vi) Braneworlds: Inflation has been implemented in 
5-dimensional phenomenological braneworld models 
lU [s^. Gravitational waves at direct detection scales 
cross the Hubble radius at high energies {Hi ^ 1, where 
£ is the bulk curvature) , hence tUg^j is directly affected by 
extra-dimensional physics. An enhancement effect in ujg^ 
arises through the modification of the Friedmann equa- 
tion, whereas a damping effect occurs via the mixing of 
massive Kaluza-Klein modes with the massless graviton 
84, 85, 86]. At direct detection scales, it is conceivable 
that these two effects cancel [stIIssI . 

Finally it is worth noting that we have ignored astro- 
physical sources, most notably from inspiralling binary 
systems of white dwarfs, neutron stars or black holes 
which could produce a significant background at frequen- 
cies of 1 mHz to 1 Hz. These sources must be subtracted 
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Phenomena 



1. Slow-roll inflation 



2. Pathological potential 



3. Bubble nucleation 



4. Turbulence 



5. Cosmic strings 



6. Pre-Big Bang / Cyclic 
models 



7. Braneworlds 



Key parameters 



Inflationary energy scale. 

'Breaking' scale(s). 

Sharp changes in V,V',V" etc. 

Bubble velocity. 
Energy scale of transition. 
Time scale of transition. 
Efficiency of energy conversion. 
+ SUSY parameters 

Characteristic turbulent scale. 
Damping scale. 
Energy scale of turbulence. 
Time scale of turbulence. 
Efficiency of energy conversion. 
+ Detail of helical/magnetic fields 

String tension. 
Average loop size. 
Intercommutation probability. 
Burst rate. 

+ Detail of loop distribution 



Detail of stringy epoch? 



Bulk curvature scale. 
Radion potential. 
+ Geometrical setup 



Ljgw(l mHz - 1 Ifz) 



< 10- 



< 10-15 



10" 



< 10-^ 
" - 10- 



< 10-^ 
10-" 



< 10"' 
10-^" 



< 10- 



References 



Eq. lUll 

[3, 50] 



55,|56 

El 



TABLE II: Summary of some possible cosmological sources of primordial gravitational wave background in the frequency range 
of future direct detection experiments (/ ~ 1 mHz - 1 Hz). Infiation may be accompanied by some (or all) of phenomena 2-5, 
while phenomena 6 and 7 are alternatives to the inflationary scenario. 



to high accuracy [s^ |^ to achieve sensitivities of 
LOgw ^ 10~^^ necessary to test inflation, and may ulti- 
mately limit direct detection experiments. The sensitiv- 
ities of the post-LISA experiments quoted in Table I de- 
pend on the usable frequency range and are significantly 
lower if frequencies 0.2 Hz are contaminated by a high 
background from unresolved white dwarfs. 

VI. CONCLUSIONS 

The generation of tensor modes is a key prediction of 
inflationary models and has yet to be confirmed by ex- 
periment. A large experimental effort is underway to 
detect a tensor mode signature in the polarization of 
the CMB. On a longer time-scale, a number of direct 
detection experiments have been proposed to detect a 
stochastic background of gravitational waves at frequen- 
cies in the range 1 mHz - 1 Hz. However, since the spatial 
scales probed by direct experiments are some 15 orders of 
magnitude smaller than the scales probed by the CMB, 
extrapolating between these scales is highly model de- 



pendent [1I3|. 

In this paper, we have used the inflationary flow equa- 
tions to assess the accuracy of extrapolating between 
CMB and direct detection scales for single-field inflation- 
ary models. Our main results are shown in Figures 1 and 
3. For models that satisfy the observational constraints 
on Us and dng/dhik, we find a conservative upper bound 
of i^gw ^ 1-6 X 10~^^. However, as shown in Figure 
most of our models have much lower values of Wgw, and 
only a small minority have ^ 3 x 10""'^^. A direct 
detection experiment with a sensitivity of ~ 10~^^ 
is therefore limited to testing a range of single field infla- 
tionary models in which the Hubble parameter is roughly 
constant between CMB scales {N « 60) and direct de- 
tection scales {N « 20), followed by an abrupt decline 
thereafter. Examples of such trajectories are shown in 
Figure 

We have also identified a number of cosmological 
sources of stochastic gravitational wave background ac- 
cessible to direct detection experiments (Table II). In 
some cases, the predicted amplitudes are far in excess 
of those generated during inflation. A high value of Wgw 
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from, say, cosmic strings produced at the end of brane 
inflation might easily overwhelm the contribution from 
tensor modes generated during inflation. In more general 
scenarios, therefore, it may be difficult for direct detec- 
tion experiments to constrain the inflationary phase even 
if experiments can achieve 'Ultimate DECIGO' sensitiv- 
ities of Wgw ^ 10~^°. 
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